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Since theoretical and computational difficulties often preclude exact 
solution of telephone trunking network problems, approximate methods 
are naturally used. A typical approach is to determine link* blocking 
probabilities and from them calculate point-to-point blocking probabilities 
by invoking independence assumptions. Although the link blocking prob- 
abilities may be quite accurate, the point-to-point blocking probability 
calculations will, in some cases, suffer from the independence assumptions. 
This paper presents a method of taking dependence into account for certain 
networks by approximating conditional probabilities which reflect the 
dependence. The approximations avoid the problems of dealing with the 
large su?ns associated ivith this problem. 

I. INTRODUCTION 

The exact analysis of telephone trunking networks often leads to severe 
computational problems due, e.g., to the large number of possible states. 
Approximate methods are thus naturally used. There has been much 
success in approximately calculating link blocking probabilities but less 
in determining point-to-point blocking probabilities. Errors in the point- 
to-point blocking probabilities can be caused by independence as- 
sumptions. 

The purpose of this paper is to take advantage of existing techniques 
for approximating link blocking probabilities (which are quite accurate) 
and develop an approach for taking link dependences into account. In 
particular, we present a method for approximating the appropriate 
conditional probabilities for cases where the traffics are Poisson (or close 
to Poisson). The extension of the approach to the case of distinctly non- 
Poisson processes, such as arise in overflows, will be reported in Ref. 1. 
To obtain point-point blocking probabilities for non-Poisson processes, 



* We use link to denote trunk group. 
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one must also solve for the blocking seen by individual traffics when 
more than one traffic is offered to a link (the equivalent random method 
gives the blocking seen by the combined traffic). This is also treated in 
Ref. 1. 

To clarify the role of independence assumptions, suppose there is 
common traffic on links 1 and 2. Let A y and A 2 be the events that the 
common traffic is blocked on links 1 and 2, respectively. Then the 
probability of being blocked on either link 1 or on link 2, P{ A x \J A 2 ) = 
P{Aj} +P{A 2 ) - P{A,C\A 2 }, is often approximated by P{ A x UA 2 | 
^P{A X ) +P[A 2 ) -P{A 1 \P{A 2 }oTevenbyP{A l VJA 2 }g*P{A 1 } + 
P{ A 2 ) . The last approximation is clearly accurate if Pf-Ax } -\- P[A 2 ] » 
P{Ai C\ A 2 ) and the first is accurate if P{A y \ A 2 ) ^ P{A 1 ) or if both 
P\Ai C\ A 2 ) and Pf-A^Pf^} are relatively small. On the other hand, 
it is easy to give examples where the neglect of dependence leads to 
non-negligible errors. 

The errors due to neglecting dependence naturally depend strongly 
on mutual traffic. To see this in a transparent case, let us reconsider 
the two-link network mentioned above with each link having the same 
number of trunks. P{Ai \ A 2 \ could vary from PJA^ (when the 
mutual traffic is zero) to unity (when there is no traffic on link 2 not 
shared with link 1). In the latter extreme case, assuming independence 
would give P{A l \ + P\A 2 ] (1 - PiA,}) for P\A Y U A 2 \ compared 
to the correct answer P{A 1 \. Thus, assuming independence could, 
in this case, conceivably overestimate the point-to-point blocking 
probabilities by something approaching 100 percent. In fact, it is 
possible to overestimate a point-to-point blocking probability on m 
tandem links by almost as much as m times by assuming independence. 

II. BASIC APPROACH 

For simplicity of explanation, first consider the situation shown in 
Fig. 1. A,, X 2 , and \ 12 are the parameters of mutually independent 
Poisson processes. Holding times are mutually independent exponential 
random variables with unity mean (or the mean is the time unit) 
here and throughout this paper. Also, throughout the paper we shall 
assume that lost calls are cleared and that the system is in equilibrium. 
Links 1 and 2 have N x and N 2 trunks, respectively. Assume iVj ^ N 2 . 

Let A i and A 2 denote the events of link 1 and link 2 being blocked, 
respectively (i.e., there are iV\ calls up on link 1 and N 2 calls up on 
link 2). 

All of the equilibrium state probabilities may actually be expressed in 
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closed form using, e.g., the form of solution given in Ref. 2, Section 7. 
For example, 

>i., Ifx-ix, X*. If, -It, »*. 

P{A,U A 2 ] = ■"- ot '" ! , ' =0 ,.'' ! , a " i ^ *' ,! , (1) 

it+u?sN x ii\ fa I i\2 ! 

where the indices are nonnegative integers (throughout the paper). 
However, the use of such exact results very quickly becomes impractical 
for all but the smallest problems even for large computers because of the 
number of computations required. 

The link blocking probabilities can usually be quite well approximated. 
For example, the following reduced load equations, 

P, = B(Nt , X, + X 12 (l - P 2 )), (2) 

P 2 = B(N 2 , X 2 + X 12 (l - P,)), (3) 

P[A l )&p lt (4) 

P[A*}S*P 2) (5) 

(B(N, X) is the Erlang B formula) are often sufficiently accurate as they 
stand.* See Ref. 3 for much more discussion of link blocking probabilities. 
Once the link blocking probabilities are determined (however they 
are determined) all that is left to calculate is P\Ax \ A 2 \ in order to 
determine P{ A, \J A 2 } . We first write down this probability exactly and 
then show a simple practical approximation. (Approximating the sums 
becomes even more important when considering more complicated 
dependences as in Section III.) To this end, note that 
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* It may easily be shown that these Pi and P 2 agree exactly with P\A,} and 
PM») to first-order terms in \«. Although the P, usually overestimate the P\Ai), 
they can underestimate (seriously, in extreme cases) when X (2 is very large. 



2650 THE BELL SYSTEM TECHNICAL JOURNAL, OCTOBER 1971 



X 



12 



X, _^_ \z_ 

Fig. 1 — Two-link network. 



where 

A12 + X 2 

q - 1 - P. (8) 

Observe that (6) is the ratio of two sums each of which is essentially 
a Bernstein polynomial (see Ref. 4). However, a Bernstein polynomial 
is of the form X^-o OpV'/^AO but our / is actually a function of i 
or of n(i/n). Each sum is also the expectation of a function of a bi- 
nomially distributed random variable. Using the following interpolation 
formula (see Ref. 5, p. 178), 

m = m + (i - oiiw + 1) - fdt - i)] 

+ (» - D 2 hW + 1) - 2/(1) + /(* - 1)], (9) 
and taking expectations yields 

g (jpT'/W) = /(»P) + ^f IK»P + 1) ~ 2/(np) + /(np - 1)]. (10) 

(10) is the usual expression for the mean and variance approximation to 
an expectation of a function of a random variable but with the second 
derivative replaced by a central difference. Using (10) on the numerator 
and denominator of (6) yields 

P{A, | A 2 \ ^P a [A, | A 2 ) 

i _l Km f Xi _ 9 j- & - ^ 1 

= "*" 2 LiVi - AT 2 y> + 1 ^ X, J 

1 , Wflg f X, "I 

£(#, - .N>, X,) + 2 |jv\ - N 2 p + 1 'J 

Interpolation, such as given on page 571 of Ref. 3, may be used to 
evaluate the Erlang B formula for a nonintegral number of trunks or an 
approximation for Erlang B in integral form can be used. 

Note that we use the same type of approximation on the numerator 
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and denominator, which are closely related, so that the ratio can be more 
accurate than either taken separately. 

(11) shows quantitatively what we expect qualitatively, namely, that 
the dependence effect gets small as the common traffic gets small. In 
particular, when X I2 = 0, P\A X \ A 2 \ = B(N U X,) = PIA^ and when 
\ 2 = 0,P\A l \A 2 ) = B(Nj - N 2) Xi) (recall that we assumed iV, ^ N 2 ). 

To examine the approximation, we find it convenient to first consider 
P a {AJ|i4 2 } where superscript c denotes complement. Actually, since 
P[A t U A 2 \ = P(A:| + P\A[ | A 2 \P\A 2 \, P\A\ \ A 2 \ is really the 
crucial quantity. With N x ^ N 2 , we have 



-1 + *" 



P a {A\ | A 2 ) = 1 -P a \A, \A 2 \ 

T m\ 1 _ Mi ~ JM l 

N 2 pq f X, 9 i 

2 IN, - N 2 p + 1 "J 

(12) 



BjN, - N 2 p, X, 

1 

B(N 1 - N 2 p, X,) 



+ 



If we let Ni = N 2 = N, then 

1 t + j^ fr.gll 

i ■ Npq f 



P a {A\\A 2 \= B(Nq ' Al) — r : l - -^f. (13) 



5(A^5, Xx) + 2 LiVe + 1 



■ 



(12) could have been derived directly using a mean and second central 
difference approximation just as (6) was approximated by (11). The 
terms multiplied by N 2 pq/2 in both the numerator and denominator 
are then the second central difference terms. If we hold Xi and p fixed 
and let N get large, it can be shown that the second central difference 
terms get small compared to the mean terms suggesting that the ap- 
proximation is accurate for large N. However, since P a {Ax \ A 2 \ — > as 
N — > oo (with X, and p fixed), investigation of this type of convergence is 
of limited practical value. 

It is probably of more interest to examine the approximation for 
fixed p and to let X, get large with N. To this end, let X : = kN with 
A- > q. For sufficiently large n and a > n, B(n, a) = 1 — n/a so that 

B{Nq, kN) =* 1 - f ■ (14) 

Hence, we obtain for large N, 
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Npq I" kN _ 2 ,Nq] 



JVgg 



-i 2 
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_#<? + 
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^ 5(ATg, fetf) (15) 

which has the interesting interpretation that with link 2 full, the number 
of trunks on link 1 to handle the Xx traffic is reduced by the average 
number of calls on link 2 which are common to link 1 (we shall elaborate 
on this below). Although (15) is intuitively appealing and obviously 
correct for p = or 1, its accuracy should be examined for p e (0, 1). 
Observe that for large N, the variance terms dominate. 

By elaborating on the interpretation alluded to above, we can see 
why B(Nq, Xi) should tend to overestimate P\A l \ A 2 \. With D, the 
event that i calls from the X 12 traffic are up, 

D { = {zXia calls up}, (16) 

observe that 

P[A t | A 2 ] = £ P\A X | A 2 C\ Di)P{Di | A 2 ] (17) 

= EBW-iiWiDi \A a ). (18) 

i = 

Now, if the X, 2 calls were never blocked on link 1, then P{Z>, | A 2 ] = 
G)pY~'- With Ni = N 2 = N, ^ = kN, it may be shown that 
XSo B(Nt - i, XO C)pY"'' -» B(N 1 - N 2 p, X.) = B(Nq, X) as N -> °o . 
B(Ni — N 2 p, XO is the blocking probability on link 1 when the number 
of trunks is reduced by the conditional mean of trunks occupied by 
X 12 calls. Hence, under these conditions, B(Nq, X,) approximates 
P\Ai | A 2 ) by ignoring the blocking of X 12 calls on link 1. If we take 
this blocking into account, we would expect that for large N, P\A l \A 2 ) 
< B(Nq, X,). This will be seen to be the case below. 
Another approximation for P{A X \ A 2 ] is 

P t \A, I A,\ - + prj)V „. /(1 _ rj) V2 ^ 2 I »> 



where 



Vp^q- Nip + N 2 pf + 4XiA r ipg - X,g + gig - # 2 p 
7-1 2X lP 



(20) 
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k - Nt - N 



pri 



= N x 



q + pr v 



(21) 
(22) 



\? i- pr 2 

and r z is the solution between and 1 of the cubic 

X,pr* + (N 2 p - p - 2V,p + X,g - X,p)rJ 

+ (-Xifl + tf.p - A^ l9 - tf 2P - g)r 2 + JV,g = 0. (23) 

This approximation is due to D. L. Jagerman (see the derivation in the 
Appendix) . 

For Ni = N 2 = N, \i = AW, and iV large, 

\/k 2 q 2 + Ak-pq — kg 
2kp 



P b {A t \A 9 )9*1 - 



(24) 



Table I shows some numerical values of the exact and approximate 
conditional probabilities along with the values of the approximations 
for large N . It is seen that, although P b \A X \ A 2 \ does not always surpass 
the simpler P a {A! \ A 2 \, it is, on the whole, superior and its behavior 
with increasing iV is clearly better. We dwelled on P a {A i | A 2 } because 
it is useful in many cases and the mean and variance approach is easily 

Table I — Numerical Values 
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extendable to cases where extreme accuracy is not required (as in 
Section III). It can be easily seen that a percentage error in the con- 
ditional probability typically causes a smaller percent error in the 
point-to-point blocking probability. 

III. MORE COMPLICATED DEPENDENCES 

The extension of the basic idea of Section II is often rather straight- 
forward. For example, consider the network in Fig. 2 which is com- 
plicated enough to illustrate the ingredients of a general approach. 

Assume (for the sake of being concrete) that AT, ^ N 3 ^ N 2 . All of 
the traffics are mutually independent Poisson. Let A„ A 2 , A a represent 
the events of links, 1, 2, 3 being blocked, respectively. Then the fol- 
lowing are approximations to conditional probabilities: 

^ V i Pl2qi2 (AT, - iV 
PIA I A \ ^ ' =0 ^ — (25) 

t =o \ * / 1-0 ]• 

v* ( N A i *,-< *?'~' 

2L, I ■ tt>32?32 (N _ ••,, 

P\A. \ A.) & >*±jfrr ^-f, (26) 

.-0 \ * / i-0 J» ! 

Z, l .' )Pl3?13 



• //^I3^13 /»r _ .XI 

p{A t i a 3 } s ':; x , » ; -Sr-ft 1 , ( 27 ) 



2^ I i JPi3?i3 2, "7T 

,=0 \ * / 1-0 J! 



where 



X12 ~l~ X 13 (l " 3 ) 1 , OQ v 

Pl2 = X 12 + Ml - P.) + X 2 + X 8 ,(l - P,) ' ? 12 - 1 " fto , (28) 



N, N 2 _ Nj 

— : o o 

_Xi_ ^ * 2 > < x ; 

Fig. 2 — Three-link network. 
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„ X 23 + X„(l - PQ _ 

P32 " x 23 + x 13 (i - p,) + x 2 + x 12 (i - p.) ' 932 - x P32 ' (29) 

„ X.a(l - P 2 ) „ _ , ™ 

Pl3 - X„(l - P 2 ) + X 3 + X 23 (l - P 2 ) ' 9i3 - 1 - fh. , (30) 

xT = x, + x 12 (i - ^AlAl^) f (31) 

and the P, are approximations to the link blocking probabilities. In 
(31), the approximation to P{A 3 \ A 2 \ is used (and we are obviously 
assuming P 3 > 0). 

The rationale behind (25)- (27) is based on (6) at least in the case 
where the P, are small so that the carried traffics are not too far from 
Poisson. We have already discussed approximating the sums of (25)- 
(27) in Section II. 

But for the network on Fig. 2 we should also calculate P{A r \J 
A 2 U A 3 ], the blocking probability seen by the X 13 traffic. Since 

PjA.U A 2 \J A 3 \ 

= P{A l } + P{A 2 \ + P{A 3 ] 

- P{A, | A 2 }P{A 2 ] - P{A X | A a }P{A a ] 

- P{A 3 | A 2 )P{A 2 ] + P{A, H A 3 | A 2 )P{A 2 ), (32) 

we see that we need only discuss the evaluation of P{A { C\ A 3 | A 2 ), 
the other quantities already having been treated. 
Just as we derived (6), we can show that (exactly) 

P(4,n A 3 I A 2 \ 

y-> iV 2 ! pV 2 *p' 2 Vp\z*p l 2 X 3 3 ~'"~' 1:| Xf'~'"~'" 

= .•„ + .-„-Mf,f .-.-jy. *i2? gggj jia\ ijl (N 3 — i 23 — i l3 )\ (JVj — i l2 — j ia )\ 

y AMp^'Pa'pia'pi' A ' 3 ~U^"' 13 X 3 ' ^-yr"» xj' 

.„+.„4f,+,,-jv, tisltaliulfe! ,£o i a ! ,£o t'i! 

(33) 

where now 

X !2 + X 23 + X13 + X 2 
X J2 + X 23 + X13 + X 2 
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A12 + A 2 3 ~T" A13 + X 2 
A12 "T A23 T A13 + X2 

These sums may be approximated by recognizing that they are both 
expectations of functions of multinominal random variables. Thus, 

E[f(i 23 , lis , *ia)] = fiha , ?13 , *ia) 

+ i[[E(in ~ W 2 *m + E(i 13 - l l3 ) 2 fi? a + E(i l2 - i 12 ) 2 8 2 12 ] 

+ 2E[(i 23 — haXiis — Ii3)h8 23 n8i 3 

+ 2E[(i 23 — i 23 )(i i2 — I\ 2 )h8 23 n8 12 

+ 2E[(im — *i 3 )(*'i2 - Ii 2 )fi8 l3 n8 12 ]}f(l 23 , i 13 , i i2 ) (38) 

where l nm here indicates expected value and where 5 2 is a second central 
difference, e.g., 

8 2 23 f(i 23 , *ia , *iz) = /(S23 + 1, lis 1 ^12) — 2/(i 23 , li 3 , in) 

+ /(S 23 - 1, t u , In) (39) 

and the operator n8 is given, e.g., by 

juWfe j *i3 » S12) = M/(*23 , *i3 + 1, H2) - f(hz , lia - 1. *!■)] (40) 
so that, e.g., 

Ii8 23 n8 13 f(l 23 , $13 j ^12) 

= M/fe + 1, *ia + 1, S12) + /(i 23 - 1, *i3 - 1, ?ll) 

- Khs - 1, Sib + 1, £12) - /(*23 + 1, lis - 1, In)]. (41) 

The variances and covariances for the multinominal distribution are 
given on page 164 of Ref. 6. Thus, 

i ik = N 2 p ik , jk = 23, 13, 12, (42) 

o) k = E(i ik - l ik ) 2 = N 2 p ik (l - p ik ), jk = 23, 13, 12, (43) 
<r%,im = E[(i ik - i ik )(iim - h»)] = -N 2 p ik p lm , 

jk, Im = 23, 13, 12(jk * Im). (44) 
One obtains 
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X3X1 



X 3 Xi 



(A' 3 + 1)(JV, + 1) 



1 2 

"T 0-13,1 



iMi(gi - 1) 



_ 9 



2 + 77T 



X3X1 



X 3 x; ' (A' 3 + 1)(AT, + 2)(A\ + 1)/. 

+ lkL.,3(5 8 - 1 (-2)57 , (-l) - 2B; 1 B; 1 + S 3 - 2 (2)57 , (l)) 

+ ^..^'(-^'(-l) - 2BVB7 1 + B,- , (l)BT 1 (l)) 
+ a 2 i3 , l2 (B;\-l)B;\-2) - 2B^B; 1 + BV(1)B7\2))]\ 
where 

N 3 = #3 — *23 — *ia , 

A', = A 7 , — 1,2 — ll3 , 

and for j = 3 or 1, 

B? = 1/Btff, , X,), 

b;\-i) = b-' - i, 
Bj\-2) = Br - 1 - jf , 



bt(U = *r + -s 



Xj 



BJ\2) = BJ 1 + 75 



A r , + 1 



+ 7^ 



X 



AT, + 1 (N, + 2)(A', + 1) 



(45) 

(46) 
(47) 

(48) 
(49) 

(50) 
(51) 

(52) 



In many cases, ?|i, U <A 2 ^ A 3 \ rnay be sufficiently well approxi- 
mated by using B(N 3 , X 3 )B(A\ , X,) for P\A l C\ A 3 \ A 2 \ since the 
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last term in (32) is typically smaller than the preceding three which 
are, in turn, typically smaller than the first three. 

IV. DISCUSSION 

We tried the method out in a number of examples and found the 
dependence analyses to improve accuracy considerably when dependence 
is indeed important and to be useful in indicating when the dependence 
effect can be safely ignored. This investigation was actually motivated 
by a study of the use of a satellite to carry telephone traffic. (This study 
will be described more fully elsewhere.) One mode of such use, called 
"variable destination," allocates up-channels (each channel consists 
of a number of trunks) from ground stations to the satellite while all 
ground stations can receive all that is transmitted downward from 
the satellite. An equivalent network for five ground stations (the outer 
nodes) is shown in Fig. 3. Between each pair of ground stations there 
is an offered traffic; thus, there are 10 (2-way) traffics. For various 
choices of the traffics and allocations, the dependence effect becomes 
significant and it was important to know when it needed to be taken 
into account and when it could safely be ignored (which leads to sim- 
plification in the study, particularly the optimization algorithms — an 
aspect of the study was to optimize the channel allocations to minimize 
the traffic blocked). The method of this paper was found quite useful 
in this regard. It should be noted that, although the N ( are very large 
and results of the form of (1) lead to huge sums, the method was com- 
putationally simple. 

The approach to dependence given here is one of approximations 
motivated by exact results in simple contexts. The results obtained 
thus far are promising and encourage further work in investigating 
the approximations and in extending the approach. Further experience 
is needed to determine the best way to apply the results to large net- 




Fig. 3 — Equivalent network for satellite variable destination mode. 
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works (see the last comment in Section III). Observe that improved 
accuracy (and, possibly, satisfactory accuracy from a practical view- 
point) in large networks can be achieved by considering only two-link 
or only two- and three-link dependences (already discussed in this paper). 
An inherent assumption that has been used is that the traffics are 
Poisson (or close to Poisson). As mentioned in Section I, the analysis of 
point-point blocking probabilities in conjunction with the use of the 
equivalent random method will be reported in Ref. 1. 
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APPENDIX 

In Section II, two approximations for P{A V \ A 2 \ were given. This 
Appendix derives the sum approximations for P b { Aj | A 2 \ . In particular, 
approximations for the following sums will be obtained:* 

The random variable £ is distributed according to the binomial law. 
One has 



x v -' 1 



-A' + j-l „X 



z 



(N-,V-2*iT' e ^ (55) 

N-i .( -. f. Xz 

Ss-si^-'r^* (56) 

in which the integration is on a circle around the origin as center. 
Let z = re ,e , then 



^2 

(N - j) 



— 1)1 Zir J -x 



(57) 

-(AW) 



ETr = r/ P'"-"'""! Tede. (58) 

J=0 t\ zt j- r i — re 



* The derivations to follow are due to D. L. Jagerman who prepared this Appendix. 
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Computing the required expectations from (57) and (58), one obtains 

j x = L f e w-™-*'* W _ i i nr) d9 (59) 

Air J--* 

2x J- T \ 1 — re I 

in which the characteristic function <£(t) is given by 

*(r) = W*\ = fa + pe iT ) n . (61) 

Thus 

Jl "a?r/ lC-)<w ' (62) 

Jm ~hi'/* Wd9t (63) 

in which 

fcx(fl) = W - i\Tt0 - 2V In r + n In (g + pre'"), (64) 

ft 2 (0) = Are' 9 - iW0 - N In r + n In (q + pre*' 9 ) - In (1 - re <e ). (65) 

The expansions of hi(0) and /i 2 (0) in powers of 6 are 

/t,(0) = \r — N hi r -\- n hi (q -\- pr) 

+ L-N + -X~\l 

\ q + pr/ 

2 



-!H»-f-4-f »'+••■■ (66) 

2 i q + pr \q + pr/ j 

/t 2 (0) = Xr - Win r + n In (g + pr) - In (1 - r) 

+ (xr - N + -^- + r -2 —)^ 
\ g + pr 1 — r/ 



(67) 

Using the method of Hayman (see Ref. 7) the functions hi(d), h 2 (d) 
will be made stationary at 6 = by defining r x , r 2 as follows: 

\ ri -N + -fp- = 0, (68) 

g + p/-i 
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Xr 2 - N + -*&*- + -r^- = 0. (69) 

q + pr 2 1 — r 2 

Thus the values of r satisfy 

Xprl + (Xg - iVp + np)r, - Nq = 0, (70) 

XpH> + (ftp — p — Np + Xg — Xp)r 2 

+ (-Xg + Np - Nq - tip - q)r 2 + Nq = 0. (71) 

For r, , the following explicit formula is obtained: 



V(Xg - N p + np) 2 + 4iVpg - \q + Np - np _ . 

r, = 2X7"^ (72) 

For r 2 , the root between zero and one must be used. Let 

ft = Xri + -SB- - »(JL)* (73) 

ft = Xr 2 + „ -p- - «(-£*-)' + r^- + (r^-)\ (74) 
g + pr 2 \g + pr 2 / 1 - r 2 \1 — r 2 / 

then 

fc,(0) = Xr, - AT In r, + n In (g + pr,) - *|M a + • • • , (75) 

h 2 (6) = Xr 2 - N In r 2 + n In (g + pr 2 ) - In (1 - r 2 ) - £/3 2 2 + • • • . 

(76) 

Since 

f [* e-*'»d6&-±= (77) 

in which the approximation becomes the more accurate the larger /3 is, 
one now obtains 



-.v, . *„ e 



Xr. 



J.^rrOz + pr,)"^^, (78) 

J 2 S **(« + pr 2 )" ^ -— • (79) 

(1 - r 2 ) V 2tt^ 2 
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